Dissipative dynamics of vortex arrays in anisotropic traps 
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We discuss the dissipative dynamics of vortex arrays in trapped Bose-condensed gases and analyze 
the lifetime of the vortices as a function of trap anisotropy and the temperature. In particular, we 
distinguish the two regimes of the dissipative dynamics, depending on the relative strength of the 
mutual friction between the vortices and the thermal component, and the friction of the thermal 
particles on the trap anisotropy. We study the effects of heating of the thermal cloud by the escaping 
vortices on the dynamics of the system. 



I. INTRODUCTION 



Recent experimental evidence [gj-^l for vortices in 
weakly interacting Bose-condensed gases stimulated a 
number of publications on structure of a rotating state 
of superfluids (see Q for a review) . Latest state of the 
art experimental technology makes possible creation and 
study of sufficiently large vortex arrays as well as 
non-destructive measurements of time evolution in a sys- 
tem with a single vortex Such possibilities allow, in 
principle, experimentally clarify the nature of interaction 
of vortices with thermal excitations. This issue plays an 
important role for understanding of vortex dynamic in 
numerous related systems, ranging from superconductors 
to superfluid cores of neutron stars (see |Q for a review) , 
and still remains a controversial issue in the literature 

In our recent papers we investigated the dissipative 
dynamics of a single vortex Q and a large vortex array 
[ fiO| in a static (non-rotating) trap at finite temperatures. 
By this we assumed a very fast relaxation of the thermal 
component towards its static configuration due to, say, 
the anisotropy of the trap. Since actual experiments are 
performed in traps, which can be made nearly perfect, the 
rotation of the thermal component can persist for a very 
long time. Moreover, in the case of a large vortex array, 
the decay of vortices reaching the border of the conden- 
sate may be accompanied by a substantial heating of the 
thermal cloud, which, in turn, affects the strength of the 
friction forces acting on the remaining vortices. Thus, 
in order to understand an experiment in realistic condi- 
tions and use its results for a fine test of vortex dynamics 
theories, one has to have a solid description of dissipative 
dynamics including various temperature effects as well as 
the role of the trap anisotropy. 

In this Letter we develop a simple, but yet a feature 
rich, model of dissipative dynamics of large vortex arrays 
in realistic anisotropic traps at finite temperatures. We 
analyze the relaxation kinetics of a system consisting of 
the rotating thermal component interacting with a con- 
densate the vortices. We show that the trap anisotropy, 
or any other mechanism of angular momentum dissipa- 
tion, is absolutely crucial for the onset of the dissipative 
dynamics. On the basis of our analysis we distinguish the 
two regimes of the evolution. The first one corresponds 



to the previously studied static case |10|, i.e. to the sit- 
uation when the trap is sufficiently anisotropic and the 
rotation of the thermal cloud quickly ceases. In this case 
the vortex array decays in a non-exponential fashion (a 
power law for the number of vortices at sufficiently long 
times). Remarkably, in this limit the lifetime of vortices 
does not depend on the parameters of anisotropy. In the 
opposite case, i.e. when the trap is nearly perfect, the 
friction of the vortices on the thermal cloud is stronger 
than the the friction of the thermal cloud on the trap 
asymmetry. This makes the vortices and the cloud stick 
to each other: the condensate rotates with the same an- 
gular velocity as the thermal particles and the angular 
velocity of the both components decreases exponentially 
with a time constant explicitly depending on the param- 
eters of the trap asymmetry. 

We discuss the heating of the thermal component due 
to the vortex array decay. The developed model, to- 
gether with certain simplifications, allows us to obtain a 
simple analytical description of dissipative dynamics and 
compare the results with the experimental data ||]. We 
believe that the presented model can be used for indirect 
study of vortex arrays dynamics based on the thermal 
component thermometry (proposed earlier in |l0|] ). 



II. THE DESCRIPTION OF THE MODEL 

We consider an evolving vortex array in an anisotropic 
trap characterized by the radial and the longitudinal fre- 
quencies LJp and LOz- The equilibrium vortex configura- 
tion in a rotating trap (in the experiment the this 
corresponds to the situation when the stirring beam is 
on) is considered in | pl| . Further on we assume that the 
initial angular velocity of the trap rotation ^Iq ^ fie, 
where J7c is the critical velocity corresponding to the ap- 
pearance of the first vortex. In this case the equilibrium 
number of vortices iV„ 3> 1, and their spatial distribution 
is uniform. 

As we will see, the relaxation of vortex arrays occurs 
on a time scale, exceeding the relaxation time in the ther- 
mal cloud. Hence, we may consider the evolution of the 
system in quasiequilibrium approximation. This means, 
that at any given moment of time, the thermal cloud 
is characterized by its equilibrium density profile, cor- 
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responding to a certain temperature T, and its angular 
velocity of the uniform macroscopic rotation Q (see Fig. 
|l|). To derive the equations of motion for the whole sys- 
tem, consisting of the condensate with the vortices and 
the thermal cloud, we first write down the conservation 
laws. Assuming that the magnetic trap does not dissi- 
pate energy, we write down the conservation of the total 
energy 



dt 



0, 



(1) 



where Ey is the energy of the condensate with the vortices 
(see below), 



E{T) - 3r^iV(r), 



(2) 



is the internal energy of the non-condensed particles, 

.C(4) 



I{T) = mR^N{T) 



C(3)' 



(3) 



is the momentum of inertia and Rt = {2T/mijjp) 
the radial size of the thermal cloud. Here 



2U/2 



IS 



(4) 



is the number of thermal particles, Co — (wzW^)^/'^, 

Tc w(iV/C(3))i/^ is the BEC transition temperature, 
and TV is the total number of particles in the trapped 
gas. Eqs. (|^),(||) imply that T :$> fi, where ^ = nog 
is the chemical potential, g — ATrti^a/m, a is the two- 
body scattering length, and m is the mass of the gas 
particles. At lower temperatures, T ^ fi^ the size of the 
thermal cloud matches the size of the condensate and, 
hence, the temperature measurements become very dif- 
ficult. Accordingly, in our quantitative discussions we 
confine ourself to the case T ^ /i. 




FIG. 1. The schematic view of a condensate with vortices 
in a rotating thermal cloud. 



In the Thomas-Fermi limit (/it ^ Tt-uOz.p) the condensate 
is characterized by the longitudinal Rz = (2/i/mti;p)^/^ 
and the transverse Rc — sizes, so that 

Rc.z 3> let where Ic = (?i^/rn/i)^/^ is the correlation (heal- 
ing) length in the condensate. The velocity of the super- 
fluid flow, generated by a rotating vortex grid, imitates 
the flow of a normal liquid with the angular velocity fi^ : 
Vs — [r2„,r]. The corresponding density of the vortex 
lines is given by |Tl| ] 



TTK ' 



(5) 



where k = h/m is the vortex circulation (vortices with a 
larger circulation are unstable (l^). The angular velocity 
fly must not be very large. Indeed, the maximum velocity 
of the superflow, which coincides with the velocity of the 
superfluid on the border of the condensate ~ ^vRc must 
not reach the velocity of Bogolyubov sound Cs — i/^/m, 
i.e. 



« if 



(6) 



since otherwise the superflow becomes unstable with re- 
spect to the spontaneous creation of excitations already 
at T = (see Landau arguments in the nature of 
the rotational friction in such a situation is an interesting 
question by itself and is not considered here, see ||l^ for 
a detailed discussion). Since the angular velocity is 
directly related with the density of the vortex array (see 
Eq. (^), the condition (^ limits the maximum number 
of vortices in a stable vortex configuration confined in a 
harmonic trap of a given frequency: 



Ny 



TTK nUJn 



»1, 



(7) 



which is a large number in the Thomas- Fermi limit. 

Since the density of the condensate energy is logarith- 
mically singular in the vicinity of the vortex line, the 
energy of the vortex array, Ey, consists of two parts: 



E,, 



d rps 



En + ly 



(8) 



where Eq is the part of the vortex energy independent on 
the velocities of the vortices, and amounts to Ny times en- 
ergy of a single vortex at rest, Eyo = nmnoRz log(-Rc/^c), 
see |ll[]. The quantity ly = mnoRzTrR'^/2 can be called 
as the momentum of inertia for the condensate with the 
vortices (moment of inertia of the vortex array). The ra- 
tio of the first to the second term in the r.h.s. of Eq.(||) 
is ~ log[Rc/lc)/Ny <C 1. This shows that in a sufficiently 
large vortex array, Ny 3> 1, its rest energy is negligible 
and the whole internal energy of the condensate coincides 
with the energy of its macroscopic rotation. 

The second equation comes from the angular momen- 
tum balance 
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-(i(T)n + AQ 



-an, 



(9) 



where My is the angular momentum of the vortex array 
and the phenomenological coefRcient a describes the in- 
teraction of the thermal particles scattered by the trap 
anisotropy. In fact, the form of r.h.s. in Eq. ^ corre- 
sponds to the dissipation of the angular momentum by 
a friction force proportional to the velocity of the gas 
cloud. This relation between the friction force and the 
gas velocity is a very general one and holds for a friction 
force acting on an obstacle in a moving media both in 
the coUisionless and hydrodynamic limits , at least for 
sufficiently small velocities. The choice of the friction co- 
efficient a depends on the details of the anisotropy model 
and requires a microscopic calculation (see below). 

Since in the reference frame rotating with the super- 
fluid the trap asymmetry moves with a subsonic veloc- 
ity (see the discussion leading to Eq.(^)), the condensate 
particles are not scattered by the trap and, therefore, no 
term proportional to the condensate angular velocity can 
appear in the r.h.s. of Eq.(^. 

Unlike in the calculation of the condensate energy (||), 
the expression for the angular momentum of a vortex has 
no singularity and, hence, the angular momentum of the 
vortex array coincides with he angular momentum of the 
superfluid 



(10) 



Eqs. (Q),(^) are still insufficient, since they comprise a 
system of two equations for the three unknown functions 
n(t), fly{t) and T{t). The formulation of the model is 
completed by the addition of the equation of motion for 
the vortices. It follows from the Magnus law with the 
friction force originating from the interaction of the vor- 
tex lines with the thermal component. In fact it is the 
continuity equation for the vortex gas moving towards 
the border of the condensate (see Appendix ^ for the 
derivation, which closely follows the arguments given in 



I): 



at 



(11) 



where 7 = 2kPsD / {k? ~\- D^), ps = mnQ and p are the 
superfluid density and the total density of the gas, and 
D is the mutual friction coefficient, characterizing the 
interaction of the vortex line with the excitations. 

Eqs.(^,(^) and ( pi] ) constitute the complete set of 
equations describing the dissipative dynamics of vortex 
arrays interacting with thermal particles. The initial con- 
dition at i = 0, when the rotation of the trap stops (i.e. 
when the stirring beam is off) , are 



r!(o) = a„(o) = r!o, r(o) = To>Ai. 



(12) 



The formulation of the model should be completed by 
ascribing a certain value to the friction coefficient a. This 



can be done, using the results of the discussion of the 
rotational properties of trapped gas, presented in [p^ . 
Consider an ideal gas, splined up in an infinitely long 
trap, characterized by the two close transverse frequen- 
cies LOx and LUy, so that lu- 



(1 ± e)^/^, where e is the 
parameter of anisotropy. Then, after the rotation of the 
trap stops, the gas continues rotation with the angular 
velocity, decreasing exponentially with the characteristic 
time constant given by the expression p[ : 



A = Re- 



4r 



'1 - 



(13) 



where Ec = l/4a;r and r is the characteristic time be- 
tween the collisions in the thermal cloud. The relaxation 
time T can be estimated as r"-"^ ~ na yjT jm, where n is 
the gas density and a is the collisional cross section. To 
generalize this result for the time r, obtained for a gas 
sample at T > Tc, for the case of a Bose condensed gas, 
we change the total density of the gas to the condensate 
density (see the argumentation in [|l5| ). 

In our model, when the condensate is absent, or the 
interaction between the vortices and the thermal com- 
ponent is negligible (this corresponds to 7 ~ I? = 0), 
the vortices and the thermal excitations evolve indepen- 
dently. From Eq.(|l|) for the temperature change in the 
thermal cloud we have 



AT 



< 1. 



(14) 



The latter inequality follows from Eq. (H) and guarantees 
that the temperature of the system does not change. In 
this approximation Eq. (H) reads 



dfl a 



(15) 



and gives the following simple decay law for the angular 
velocity of the thermal cloud 



(16) 



where ta = //a is the relaxation time related to the fric- 
tion of the moving gas and the trap anisotropy. Introduc- 
ing the requirement, that the presented model recovers 
the result of Guery-Odelin iQ, we set = A and find, 
that a — IX. 

The appearance of imaginary part in the relaxation 
rate ( |l^ ) for lot ^ implies an oscillatory behavior 
in the onset of equilibrium of a free gas cloud (see jlj] 
for the discussion). As we shall see, the evolution of the 
combined system of vortices and thermal excitations ex- 
plicitly depends on the relaxation rate (13) only in the 
so called "weak anisotropy" limit, when both the vortex 
array and the non-condensed gas rotate with the same 
angular velocity. Since, according to Eq. (^, the num- 
ber of vortices in the vortex array is proportional to the 
angular velocity fly, such oscillations would imply either 
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periodic oscillations of the vortex density (if the oscilla- 
tions of ily are sufficiently small) , or a reversible destruc- 
tion and nucleation of vortices close to the border of the 
condensate. Well below Tc, the nucleation of vortices fl^ j 
is slow compared with kinetic times (p^, and, thus, the 
number of vortices, as well as the angular velocity of the 
condensate rotation, can only decrease in the course of 
dissipative dynamics. Although the detailed description 
requires a more thorough discussion, we neglect the os- 
cillations and keep a purely real in the remaining discus- 
sion. Collective oscillations of vortices and the thermal 
cloud are reviewed in [|l6|. 



III. ANALYSIS OF THE DISSIPATIVE 
DYNAMICS. 

As discussed in the previous section, Eqs. (|l|),(||) and 
( pi] ) form a close set of equations describing the evolution 
of a large [Ny ^ 1) vortex array emerged in a rotating 
thermal cloud. Although the equations look simple, their 
analytical solution requires certain simplifications. 

First of all we qualitatively analyze the heating of the 
thermal cloud in the course of dissipation of rotational 
energy of the gas. According to Eq. (^4|) the rotational 
energy of the thermal cloud is always less than its inter- 
nal energy and therefore makes no contribution to the 
heating of the system. Thus we can neglect it in the 
energy conservation law and express the temperature of 
the thermal cloud from Eq. (^ as a function of angular 
velocity of the condensate 



where the dimensionless parameter 



1/4 



(17) 



where A is the ratio of the rotational energy of the con- 
densate to the internal energy of the thermal cloud. It 
can be found with the help of Eqs. (0) and (HG) 



A 



2E 



fig 
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ujI (noa3)i/2 



(18) 



and can well be large. The dimensionless parameter A 
indicates how important are the heating effects in the 
course of the dissipative dynamics. On the other hand, it 
can be shown that the temperature increase T/— To Tc. 
This allows us to neglect the temperature dependence of 
Ps and ly. 

Remarkably, in Eq. (^ we can neglect also the time de- 
pendence of the momentum of inertia of a thermal cloud. 
Indeed, the term containing the time derivaty of the tem- 
perature, /f2, turns to be always small compared to ly^v- 
ifl/lih, (flM/^Dii^p/^o)^ < 1- Using this fact in 
Eqs. (^) and (|lj) we obtain 



V7ri^ - rj„) + f2 

"0 



(19) 



V 



{I + Iv)jilo _ Tr 

a T„ ' 



characterizes the relative strength of the mutual interac- 
tion between the vortices and the thermal particles and 
the interaction of the thermal cloud and the trap asym- 
metry. Here = 'jVIq and r^^ = a/{I + 1^) are the 
characteristic relaxation times, determined respectively 
by the interaction of thermal component with the vor- 
tices and with the trap anisotropy. 

The parameter r] distinguishes two limiting regimes of 
relaxation. First we turn to the static trap limit, when 
Tj^. Ty (77 ^ 1). In this case at the initial stage the fric- 
tion of the thermal cloud on the trap anisotropy is much 
stronger than the friction on the condensate. Therefore, 
at this stage the relaxation of the thermal cloud should 
proceed in the same way as in the absence of vortices, 
i.e. the angular velocity f2 exponentially drops on the 
time scale ta ^ t^j <C t« , while the angular velocity of a 
condensate remains almost unchanged, Cly « 0. Neglect- 
ing the first term in the brackets of r.h.s of Eq. ( pj| ) and 
using Eq. (^ we reduce Eq. (|l^) to Eq. ( [l5| ) . This quali- 
tative picture remains valid until the thermal clod nearly 
stops and its angular velocity reaches the very small value 
< rjfly/flo ^ Hq. At this point the trap becomes effec- 
tively static and, setting 51 « in Eq. (pi]), we find the 
equation of motion for the condensate angular velocity 



(20) 



If the time dependence of the coefficient 7 can be ne- 
glected, that either for A ^ 1, or for very small times 
when Sl„ « r2o, the solution of this equation is given by 
the expression 



ny{t) 



no 



1 + t/Ty ' 



(21) 



which recovers the result previously obtained in |lO| . 

To analyze the solutions of Eq. (20) in the strong heat- 
ing limit, A 3> 1, we note that 7 cx T and rewrite Eq. 
( pO| ) with the help of Eq. (|l^) and the dimensionless 
units for the condensate angular velocity fly = fly/flo 
and the time t = </Tfl(To): 



^ = -(i + A{i-nl)y/'nl 

at 



(22) 



Integrating it we find the exact implicit solution Qy{t): 

il + A-Anlf/\ AQy ( A^l \ 

{l + A)VLy ^2(l + A)5/4^-i- 1^1 + ylj 

where Fa.b,c{x) is the hypergeometric function and the 
integration constant C is to be determined from the ini- 
tial condition a t = 0. In the considered limit of large A 
we find, that 
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c = 



v^r(3/4) 
4r(5/4)Ai/4- 



Then, for small fit,, which corresponds to large times, 
when the rotation of the condensate is almost stopped, 
the hypergeometric function can be expanded in series 
of its small argument and in the lowest order in 1/A 
we obtain the asymptotic expression for the condensate 
angular velocity 



0.4 + Vr.(r/)' 



which is valid as soon as t ^ Ty{Tf). This means, that 
the quantity Ty{Tf ) plays the role of the vortex array life 
time and is mainly determined by the friction coefficients 
at the final temperature. This fact is fairly intuitive, 
since the strength of the interaction of the vortices with 
the thermal cloud quickly increases with the rise of the 
temperature and thus most of the relaxation of the sys- 
tem occurs when its temperature is close to its final value. 
Remarkably, the dynamics of the vortex array in a static 
trap limit, that is when anisotropy is strong, turns out to 
be independent from the parameters of the anisotropy. 

In the opposite, rotating trap limit, when the 
anisotropy is weak and r^j ^ r^, (77 3> 1), the friction 
between the vortices and the thermal cloud is stronger 
than the friction of the thermal excitation on the trap 
asymmetry. As the result, the both components stick to 
each other due to their mutual interaction and dissipate 
the angular momentum together, as a single rigid body. 
Neglecting the last term in the brackets of r.h.s. of Eq. 
( |l9[ ) for the initial conditions given by (^2|) we find the 
solution f2 = Using this relation between the angular 
velocities in Eq.(g) we arrive at 



dt ~ I + Iv 



^7; 



(23) 



Again, when the heating of the gas sample is small, the 
coefficient a is time independent, and the solution of this 
equation is simple: 



(24) 



i.e. in this, "rotating trap" approximation, the power law 
( pT| ) changes into the simple exponential decrease of the 
angular velocities, with the same time constant for both 
the vortex array and the thermal cloud. 

To discuss the solutions of Eq. ( p3| ) in the strong heat- 
ing limit A ^ 1 we have to make certain assumptions 
about the trap anisotropy. For example, according to 
Eq. (|l3|) , when the trap asymmetry is small, AeujT <^ 1, 
we have a oc T^. Then Eq.(E3) has the solution 



2cxpi-t/Tn{Tf)) 
l + expi~2t/TR{Tf))- 



(25) 



This means that again, the relaxation time at the finite 
temperature TuiTf) plays the role of the vortex array life 



time. At larger times, t » Tv{Tf ), the angular velocity 
is small and the solution has an approximate asymptotic 
form 

ny = 2^0 eM-t/MTf)) ■ 

Remarkably, in logarithmical approximation the solution 
(p5| ) holds also for a situation, when the trap asymmetry 
is strong, ieujT > 1, and a cx T*^ (see Eq.(|3|)). 

Using the presented above analysis we now can draw 
a qualitative picture of dissipative dynamics. To do this, 
we introduce the dimensionless ratio of the momenta 
of inertia of the thermal cloud and the vortex array: 

= I/Iy - (To//i)4(noa3)i/2. Since in a weakly interact- 
ing gas the gaseous parameter is small: (noa^Y^^ <^ 1, 
at T fi the condensate is heavier than the thermal 
cloud, and both the rotational energy and the angular 
momentum are dominated by the condensate contribu- 
tions. Thus, in the rotating trap limit (77 <^ 1) the re- 
laxation time T/f is mainly determined by ly and is much 
larger than the spinning down time of a free thermal 
cloud TA In the static trap limit the role of the 

parameter ^ is more subtle. At sufficiently large times, 
t }^ TA, we can neglect the term itl in Eq. (^, which 
thus transforms into lyfl = —aQ. Comparing this equa- 



tion with Eq. (20) we find that Q = TyO^/fioj i-e. the 



thermal cloud angular velocity decays as a power law and 
closely follows the relaxation of the vortex array. In the 
same time the heating parameter A is related to the pa- 
rameter C: strong heating of the sample is only possible 
when ^ ^ 1. As the temperature grows, the number of 
the condensed particles decreases and at some point the 
condensate may become so small, that C ^ 1. One can 
see, that the crossover between the two regimes occurs 
at the temperature below T^i so the condensate does not 
disappear. Later on, the heating of the system stops and 
the system relaxes according to Eqs. (|2|) and (^) for 
ry ^ 1 and 77 ^ 1 respectively. 

The presented analysis shows that the anisotropy, or 
any other mechanism of angular momentum dissipation, 
is absolutely crucial for the onset of dissipative dynamics. 
It is seen already from the fact, that when th — > 00, 
Eq.(p4[) gives ^lv{t) = r2(i) = fioj i-e. the rotation of the 
vortex grid persists forever. 

Although the presented analysis is sufficient for a qual- 
itative understanding of the dissipative dynamics in the 
whole parameters range, the quantitative description re- 
quires a numerical solution of the initial Eqs. (0),(|^) and 
( ll|) . For the sake of completeness we provide the results 
of such calculation based on the experimental parameters 
of |^,||. We use the following values for the model pa- 
rameters: the condensate density uq — 2x 10^^cm~^, the 
trap frequencies ujp — 2Trx 169s~^ and uj^ = 27rx 11.7s~^, 
and the temperature Tq ~ /i = 80nK . Then, the conden- 
sate size is Rc — Aiim and, at the reported frequency of 
the stirring ilo = 27r x 135s~^, the vortex array consists 
of Ny ~ 10 vortices. 
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FIG. 2. Solid (dashed) line represents Qv (Q) obtained by 
numerical solution of Eqs. and ( [ll| ) for the param- 

eters of the ENS experiment. Dotted and dash-dotted lines 
represent approximate solution (^sj) with tr calculated for the 
strong (e ^ Ec) and the weak (e <^ ec) anisotropy limits. 

We provide the results of the numerical calculations 
for the trap anisotropy parameter e = 0.06 (see Figs. ^ 
and 1^). Our model system turns out to be well in the 
rotating trap limit (rj ~ 10) and, in the full agreement 
with our qualitative analysis, « r2„ and either of the 
angular velocities relaxes exponentially (see the solid and 
the dashed lines on Fig. ||). The solid line on Fig. ^ 
illustrates the heating of the system in the course of the 
dissipative dynamics. 
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FIG. 3. Solid line represents temperature of the thermal 
cloud obtained by numerical solution of Eqs. (§),(^, and ( |ll| ) 
corresponding to Fig. ^ Dotted and dash-dotted lines repre- 
sent temperature obtained from Eq. ( [l7| ) using approximate 
solutions presented in Fig. |^. 

The comparison of the exact numerical solution with 
our approximate analytical results is a bit difficult, since 
this value of the anisotropy parameter is close to ec (see 
Eq.(p^). The dashed and the dotted dashed lines on 



Figs. Id and ^ represents the analytical results based on 
Eq. ( Pq ) with calculated in the strong e 3> Ec and the 
weak e ^ Ec anisotropy limits in Eq.(^3|) respectively. 
We note, that all of the presented results are in a good 
agreement with the experimentally found vortex array 
life time ~ Is. 



IV. CONCLUSION 

In this Letter we formulate a realistic model, describing 
the evolution of a large vortex array in a rotating thermal 
cloud interacting with the trap anisotropy. The model 
allows both the simple analytical description of the most 
important limiting cases as well as a very straightforward 
numerical analysis. 

First we show that the angular momentum dissipation 
mechanism, such as the trap anisotropy is absolutely cru- 
cial for the onset of dissipative dynamics of vortices in 
trapped Bose-condensed gases. We demonstrate, that if 
the total angular momentum of the system is conserved, 
the rotation of the vortex array persists forever and, con- 
sequently, the number of the vortices in the condensate 
remains constant. 

Since the flow of the condensate is a superflow, the con- 
densate itself can not interact with the trap anisotropy. 
On the contrary, the thermal excitations are scattered 
by the trap and can dissipate the angular momentum. 
Thus, it is just the mutual friction between the vortices 
and the gas particles that allows the condensate to trans- 
fer its angular momentum to the excitations and, subse- 
quently, to the trap. Depending on the relative strength 
of interaction of the thermal gas with vortices and with 
the trap anisotropy we distinguish the two regimes of the 
dissipative dynamics. 

First we identify the limit, when the trap anisotropy is 
very strong and the rotation of the thermal cloud quickly 
stops. Then the vortex array interacts with the thermal 
cloud at rest and subsequently looses the vortices due to 
the mutual friction between the vortices and the excita- 
tions (static trap limit). This case has been previously 
studied in Q for a single vortex and in for a large 
vortex array. 

In the opposite limit, when the trap is nearly perfect, 
the thermal cloud stops slowly. In this approximation 
the mutual friction between the vortices and the ther- 
mal cloud makes the two components stick to each other, 
forming a single rotating object with the momentum of 
inertia equal to the sum of the momenta of inertia of the 
condensate and the thermal cloud. This composite body 
is usually much heavier than the thermal cloud and hence 
its relaxation proceeds on a time scale much larger than 
that in the Guery-Odelin scenario jl4j, where the calcu- 
lation was done above Tc and only the thermal cloud was 
taken into account. 

Another feature of a large vortex array evolution is the 
heating of the thermal component caused by the energy 
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released by the vortices escaping the condensate. We 
show, that though the heating itself can be very sub- 
stantial, it does not change the qualitative picture of the 
dissipative dynamics. Since the friction coefficients a and 
7 both grow very quickly with the temperature increase, 
most of the relaxation processes occurs when the tem- 
perature is already high, i.e. when the temperature is 
close to Tf. Hence, for estimation purposes, one can use 
a much simpler calculation with T = Const with all the 
model parameters taken at the final temperature. 
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APPENDIX A: 



DYNAMICS OF VORTEX 
ARRAYS. 



To derive the equation of motion for the vortices inside 
a rotating thermal cloud, we closely follow the derivation 
presented in |]l0| . The configuration of the superfluid is 
described by the superfluid velocity Vs and by the two 
continuous functions: the average vortex density ny{t, r) 
and velocity Vi,(i, r). Since the number of vortex lines 
is a locally conserved quantity, the introduced functions 
satisfy the continuity equation: 



dr 



= 0. 



(Al) 



The velocity of the superfluid satisfies Maxwell-like 
equations 



rotVs = 2TiK,nyZ, divv^ = 0, 



(A2) 



where z is a unit vector along the axis of the cylinder [ plj . 
These equations are solved with the boundary condition 
that the normal component of vanishes on border of 
the condensate. The vortex velocity is related to the 
superfluid velocity through the Magnus law O: 



Ps [Vt, 



(A3) 



where ps is the superfluid (mass) density and F is the 
friction force acting on the vortex line (per unit length). 

At finite temperatures, the force F originates from the 
scattering of thermal excitations from the moving vor- 
tices and consists of two terms: 



F = -D[wy - v„) + D'[kz, (v^, - v„)], 



(A4) 



where v„ is the velocity of the thermal gas at the position 
of the vortex line. Here D and D' are the longitudinal 



and transverse friction coefficients respectively. Both of 
them are temperature dependent and D' = p„, where pn 
is the density of the normal component 0. The friction 
coefficient D is also temperature dependent and ranges 
from D (X Knt'l''''T^ j [iJ I'^T? for very small temperatures 
T < ^ lei to D w 0.lKTO5/2y^i/2/^3 |] for ^ < T < 



As discussed in |10 , the vortex density is spatially uni- 
form: n^(t, r) = riy and the superfluid velo city inside the 
vortex array can be found from Eq.(A2) by using the 
Stokes theorem. It has no radial component (i.e. it sat- 
isfies the necessary boundary condition) and is given by 



V5(r) = 7r[K, r]n^ 



(A5) 



The solution of Eqs.(A3) and (A4), with from Eq.(A5) 
and v„ = [ri,r], can be represented in the form: v^, = 
^('■)f _(_ yi'i') J/?;^ f ] ^ where 



and 



D2 



(A6) 



(A7) 



are the radial and tangential components of the velo city 
field, and p = -I- p„ is the total density. Eqs. (A6) 
and (A7) show, that if the condensate rotates together 
with the normal component, Vs — w„, then v"^^ — w„ 
and the radial vortex velocity is absent, i.e. the vortex 
and the thermal particles move together and are in the 
equilibrium in the reference frame defined by the rotating 
cloud. 



Using Eqs.(A5) and 
form 



), we rewrite Eq.(Al) in the 



dn^, (t) 
dt 



where 7 = 2kPsD/{k^p^ + D'^) w 2D/ up. Finally, ex- 
pressing the vortex density through the angular velocity 
of the condensate rotation (see Eq.(||)), we obtain the 
desired equation of motion for the vortices (E^. 
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